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•  Discharge  characteristics  of  a  multicell  lithium-ion  battery  with  non-uniform  cells  are  studied. 

•  The  relationship  between  discharge  capacity  of  a  multiple  battery  and  its  topology  is  unraveled. 

•  A  mathematical  model  for  evaluating  the  discharge  capacity  of  multicell  batteries  with  general  topology  is  derived. 
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This  paper  addresses  the  discharge  characteristics  of  a  multicell  lithium-ion  battery  with  nonuniform 
cells.  By  using  a  multicell  battery  with  two  different  types  of  cell,  we  experimentally  investigate  the 
relationship  between  discharge  capacity  and  battery  topology.  In  the  series  connection  of  nonuniform 
cells,  the  total  capacity  of  the  battery  is  governed  by  a  cell  with  the  smallest  capacity.  In  the  parallel 
connection,  the  total  capacity  becomes  the  sum  of  cells  in  parallel.  In  a  multicell  battery  with  the  series 
-parallel  connection,  the  discharge  capacity  depends  on  its  configurations  such  as  connection  and 
number  of  nonuniform  cells.  Also,  we  derive  a  mathematical  model  for  discharge  characteristics  of  the 
multicell  battery  based  on  the  equivalent  circuit  model  of  an  individual  cell.  For  a  multicell  battery  with  a 
general  configuration,  the  discharge  capacity  is  successfully  evaluated  by  numerical  simulations  of  the 
mathematical  model. 

©  2013  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Multicell  lithium-ion  battery  is  a  system  of  energy  storage  that 
consists  of  multiple  lithium-ion  cells  connected  in  series  and  par¬ 
allel.  Due  to  its  high  energy  capacity,  it  is  used  in  various  applica¬ 
tions  such  as  electric  vehicles  and  embedded  systems  [1-4].  A 
multicell  battery  is  initially  synthesized  using  cells  with  a  uniform 
characteristic  [3,4].  However,  due  to  the  repetitive  charge/ 
discharge  operations,  the  initial  uniformity  of  cell  characteristics 
gradually  breaks  down.  It  is  widely  recognized  that  such  nonuni¬ 
form  characteristics  deteriorate  the  performance  of  the  battery  and 
cause  a  problem  on  safety  and  reliability  [2,4,5].  In  order  to  avoid 
the  nonuniform  characteristics,  understanding  the  influence  of 
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nonuniform  characteristics  on  the  performance  is  of  basic  impor¬ 
tance.  This  leads  to  design  and  management  of  high-confidence 
multicell  batteries  [6-9]. 

Recently,  researchers  have  focused  on  the  discharge  character¬ 
istics  of  multicell  batteries  with  nonuniform  cells.  An  effect  of  low- 
capacity  cell  on  the  battery  performance  was  mentioned  for  the 
design  of  management  systems  [4].  A  circuit-based  model  was 
derived  to  estimate  the  distribution  of  discharge  currents  in  a 
general  multicell  battery  [10,11].  An  iterative  calculation-based 
method  was  also  proposed  for  simulating  discharge  behaviors  [12]. 

The  purpose  of  this  paper  is  to  study  the  discharge  character¬ 
istics  of  a  multicell  lithium-ion  battery  with  nonuniform  cells.  As 
shown  in  Ref.  [11],  the  discharge  characteristics  become  compli¬ 
cated  due  to  current  effect,  recovery  effect,  temperature  effect, 
aging  effect,  and  their  interactions.  To  the  best  of  our  knowledge,  no 
experimental  analysis  of  the  discharge  characteristics  in  terms  of 
nonuniform  cells  has  been  reported.  Here,  from  a  practical  view¬ 
point,  it  is  important  to  clarify  the  relationship  between  discharge 
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characteristics  and  battery  configurations  such  as  connection  and 
number  of  nonuniform  cells.  The  multicell  battery  studied  in  this 
paper  includes  two  different  types  of  cell:  one  has  large  internal 
resistance  and  small  capacity,  and  the  other  has  small  internal 
resistance  and  large  capacity.  By  using  a  battery  with  one  former 
cell  and  many  latter  cells,  we  experimentally  investigate  the  rela¬ 
tionship  between  discharge  capacity  and  battery  topology.  The 
topology  implies  three  distinguishable  connections:  series,  parallel, 
and  series-parallel  connections.  Also,  we  derive  a  mathematical 
model  for  evaluating  the  discharge  capacity  of  a  multicell  battery. 
The  effectiveness  of  the  derived  model  is  demonstrated  with  ex¬ 
periments  and  numerical  simulations  of  a  battery  with  a  general 
configuration. 

The  rest  of  this  paper  is  organized  as  follows.  Section  2  describes 
basic  concepts  of  multicell  battery.  The  experimental  setup  and 
procedure  used  in  this  paper  are  described  in  Section  3.  Section  4 
experimentally  studies  the  discharge  characteristics  of  a  multicell 
battery  in  the  series,  parallel,  and  series-parallel  connections. 
Section  5  analyzes  a  dependence  of  the  battery  configurations  on 
the  discharge  capacity  in  the  series-parallel  connection  by  using 
experiments  and  numerical  simulations.  Section  6  is  the  conclusion 
of  this  paper  with  a  brief  summary  and  remarks. 

2.  Basic  concepts  of  multicell  battery 

Here,  basic  information  about  connections  and  characteristics  of 
existing  multicell  batteries  is  introduced.  There  exist  three  distin¬ 
guishable  connections  in  a  multicell  battery,  which  we  refer  to  as  (i) 
series,  (ii)  parallel,  and  (iii)  series-parallel  connections  [3,4].  Se¬ 
ries-parallel  connection  is  further  classified  according  to  the  ex¬ 
istence  of  cross  coupling:  see  Section  4.3  and  Fig.  6.  In  this  paper,  Ns 
denotes  the  integer  number  of  cells  connected  in  series  and  NP  the 
integer  number  of  cells  in  parallel.  Also,  (Ns,NP)  denotes  the  com¬ 
bination  of  number  of  cells  in  series  and  in  parallel.  A  commercial 
multicell  battery  possesses  a  mixture  of  the  three  connections. 
From  graph  theory  [13],  every  series-parallel  network  is  aggre¬ 
gated  as  either  series  or  parallel  network.  This  suggests  that  arbi¬ 
trary  connection  is  aggregated  as  either  basic  series  connection  or 
parallel  connection  [11,12].  On  the  basis  of  this  fact,  in  this  paper  we 
define  the  basic  series  and  parallel  connections  as  (Ns,NP)  =  (2,1) 
and  (NS,NP)  =  (1,2).  Also,  a  multicell  battery  is  generally  synthesized 
using  cells  with  a  uniform  characteristic.  In  this  uniform  case,  the 
discharge  characteristics  on  the  series-parallel  connection  are 
explained  in  terms  of  the  basic  series  and  parallel  connections 
[11,12].  On  the  other  hand,  in  the  case  of  nonuniform  cells,  no 
explanation  has  been  reported  in  terms  of  basic  connections. 

3.  Experimental  system  and  procedure 

This  paper  studies  the  multicell  battery  in  which  only  one  cell  is 
different  from  the  others.  In  Section  4,  based  on  the  idea  of  ag¬ 
gregation  mentioned  above,  we  experimentally  analyze  the 
discharge  characteristics  of  the  series  and  parallel  connections  with 
two  cells  shown  in  Fig.  1.  Then,  the  discharge  characteristics  are 
analyzed  in  the  case  of  series-parallel  connection  denoted  as 
(NsJVP)  =  (2,2). 

The  experimental  system  and  procedure  used  in  this  paper  are 
described  below.  Fig.  2  shows  the  overview  of  the  experimental 
system.  In  order  to  uniform  the  initial  conditions,  individual  cells 
were  first  charged  to  their  full  capacities  under  the  conditions  of 
constant  current  and  constant  voltage.  Every  cell  was  charged  until 
its  open-circuit  voltage  reached  4.15  V  ±  10  mV.  Then,  the  multicell 
battery  with  the  full-charged  cells  was  connected  to  the  electronic 
load  (TEXIO,  LW151-151D)  and  was  discharged  until  the  terminal 
voltage  of  a  cell  reached  its  cut-off  voltage  of  2.5  V  ±  10  mV.  During 


Fig.  1.  Basic  series  and  parallel  connections  of  multicell  battery  with  two  different 
cells,  (a)  Series  connection,  (b)  parallel  connection. 

the  discharge  operation,  the  electronic  load  was  operated  under  the 
constant  resistance  at  4.8  Q  per  a  cell.  The  terminal  voltage  and 
current  of  every  cell  were  measured  and  recorded  using  the  data 
loggers  (KEYENCE,  NR-500).  The  current  was  measured  with  the 
current  sensing  resistors  at  5  mQ  (Welwyn  Components,  OAR-3). 
Then,  the  discharge  capacity  was  calculated  by  numerically  inte¬ 
grating  the  measured  current  during  the  period  of  discharge 
operation.  All  experiments  in  this  paper  were  conducted  at  room 
temperature.  Ambient  temperature  for  the  experiments  is  pre¬ 
sented  in  each  caption  of  figures. 

Two  types  of  cylindrical  cells,  which  were  ICR  (BYD,  ICR18650) 
and  NCR  (Panasonic,  NCR18650),  were  used  in  the  experiments.  The 
ICR  cell  was  used  as  the  one  different  cell.  The  rated  characteristics 
of  ICR  and  NCR  are  summarized  in  Table  1.  The  nominal  discharge 
capacity  of  ICR  is  2350  mAh  and  is  smaller  than  the  nominal  ca¬ 
pacity  (2900  mAh)  of  NCR.  Fig.  3  shows  the  measured  discharge 
voltage  and  internal  resistance  for  ICR  and  NCR.  The  internal 
resistance  was  measured  using  the  i-R  drop  between  terminal 
voltage  and  open-circuit  voltage  [3].  In  Fig.  3,  the  terminal  voltage 


Fig.  2.  Experimental  system  of  multicell  lithium-ion  battery. 
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Table  1 

Rated  characteristics  of  two  tested  cylindrical  cells. 


ICR18650 

NCR18650 

Nominal  capacity 

2350  mAh 

2900  mAh 

Nominal  voltage 

3.7  V 

3.6  V 

Cut-off  voltage 

2.5  V 

2.5  V 

of  ICR  is  lower  than  NCR,  because  of  the  large  internal  resistance 
and  the  low  open-circuit  voltage  of  ICR. 

4.  Experiments 

This  section  experimentally  investigates  the  discharge  charac¬ 
teristics  in  the  basic  connections  of  nonuniform  cells.  Here,  the 
discharge  characteristics  in  the  series-parallel  connection  are 
evaluated  in  terms  of  the  basic  series  and  parallel  connections.  In 
this  paper,  A  +  B  denotes  the  series  connection  of  cells  A  and  B,  and 
A  ||  B  the  parallel  connection  of  cells  A  and  B,  respectively.  Also, 
Qs(A)  denotes  the  discharge  capacity  of  cell  A  in  the  case  of  a  single 
operation. 


Table  2 

Discharge  capacities  of  two  cells  in  series  connection  at  ambient 
temperature  of  28.3  °C.  The  symbol  f  denotes  the  value  of 
discharge  capacity  under  the  single  operation. 


Discharge  capacity/mAh 

11 

2307  (23261) 

Nl 

2299  (2839f) 

11  +  Nl 

2303 

of  2.4  Q(=4.8  Q/2).  The  results  are  shown  in  Table  3  and  Fig.  5(a).  In 
contrast  to  the  series  connection,  the  terminal  voltage  of  each  cell  is 
common  as  shown  in  Fig.  5.  During  the  discharge  operation,  each 
cell  discharges  from  an  initial  voltage  to  its  cut-off  voltage.  As  a 
result,  the  discharge  capacities  of  II  and  N2  almost  correspond  to 
Qs(H)  and  Qs(N2),  respectively:  see  Table  3.  Thus,  in  the  parallel 
connection  of  two  cells  with  nonuniform  characteristics,  the  total 
capacity  of  the  battery  becomes  the  sum  of  individual  cells.  In  this 
connection,  the  discharge  current  of  the  battery  is  not  equally 
distributed  to  each  cell  as  shown  in  Fig.  5(b). 

4.3.  Series-parallel  connection 


4.1.  Series  connection 

First,  we  experimentally  investigate  the  discharge  characteris¬ 
tics  of  two  cells  in  series.  As  shown  in  Fig.  1(a),  one  ICR  (II, 
2326  mAh)  and  one  NCR  (Nl,  2839  mAh)  were  connected  in  series 
and  discharged  to  a  load  of  9.6  Q(=4.8  Q  x  2).  Table  2  summarizes 
the  measured  discharge  capacity,  and  Fig.  4  shows  the  measured 
discharge  voltage  and  current  of  each  cell.  Discharge  capacity, 
voltage,  and  current  are  also  shown  in  the  case  of  a  single  operation. 
In  the  series  connection,  the  current  of  each  cell  is  common  as 
shown  in  Fig.  4(b).  As  a  result,  the  terminal  voltage  of  II  reaches  its 
cut-off  voltage  earlier  than  Nl:  see  Fig.  4(a).  Here,  the  discharge 
capacity  of  the  battery  is  2303  mAh  and  is  close  to  Qs(Il).  In  this 
case,  Nl  has  nonzero  amount  of  charge  that  corresponds  to 
Qs(Nl)  -  Qs(H)  at  the  end  of  discharge  operation.  Thus,  the  total 
capacity  of  the  battery  is  determined  by  the  cell  with  the  small 
capacity  in  the  series  connection. 

4.2.  Parallel  connection 

Next,  the  discharge  characteristics  of  two  cells  in  parallel  are 
experimentally  investigated.  II  and  another  NCR,  denoted  as  N2 
with  2820  mAh,  were  connected  in  parallel  and  discharged  to  a  load 


600 

450 

300 

150 
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Fig.  3.  Discharge  voltage  and  internal  resistance  of  ICR  and  NCR  at  ambient  temper¬ 
ature  of  25.0  °C. 


Here,  the  discharge  characteristics  in  the  series-parallel 
connection  are  investigated.  The  connections  of  four  cells  are 
shown  in  Fig.  6.  As  mentioned  in  Section  2,  the  connections  are 
distinguished  according  to  the  existence  of  cross  coupling.  II  and 
three  NCR  cells  (Nl,  N2,  N3)  were  connected  and  discharged  to  a 
load  of  4.8  Q. 


Fig.  4.  Discharge  voltage  and  current  of  II  and  Nl  in  series  connection  at  ambient 
temperature  of  28.3  °C.  The  dashed  lines  show  the  discharge  characteristics  in  the  case 
of  single  operation,  (a)  Discharge  voltage,  (b)  discharge  current. 
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Table  3 

Discharge  capacities  of  two  cells  in  parallel  connection  at 
ambient  temperature  of  28.5  °C. 


a 


I 


739 


Discharge  capacity/mAh 

11 

2301  (2326f) 

N2 

2833  (28201) 

11  ||  N1 

5134 

4.3 A.  With  cross  coupling 

First,  we  present  the  experimental  results  for  multicell  battery 
with  cross  coupling.  Table  4  summarizes  the  measured  discharge 
capacity  and  Fig.  7  shows  the  measured  discharge  voltage  and 
current  of  each  cell.  In  Table  4,  the  discharge  capacities  of  II  and  N2 
correspond  to  Qs(H)  and  Qs(N2),  respectively.  The  discharge  ca¬ 
pacities  of  N1  and  N3  are  around  2550  mAh.  The  capacities  are 
smaller  than  Qs(Nl)(=Qs(N3))  and  almost  half  of  the  sum 
Qs(Il)  +  Qs(N2).  The  total  discharge  capacity  of  the  battery  is 
5121  mAh.  These  results  imply  that  the  total  discharge  capacity  is 
restricted  by  II  ||  N2.  This  is  consistent  with  the  results  on  the  series 
connection  in  Section  4.1.  The  connection  of  cells  in  Fig.  6(a)  is 
aggregated  as  the  series  connection  of  II  ||  N2  and  N1  ||  N3.  From 
the  results  on  the  series  connection,  the  discharge  capacity  of  the 
battery  is  restricted  by  II  ||  N2.  Then,  from  the  results  on  the  par¬ 
allel  connection,  the  discharge  capacity  is  calculated  as  the  sum 
Qs(Il)  +  Qs(N2).  The  sum  is  almost  equal  to  the  experimentally 
estimated  capacity.  The  distribution  of  currents  in  II  and  N2  on 
Fig.  7(b)  is  almost  the  same  as  in  the  case  of  the  parallel  connection 
of  II  and  N2  in  Fig.  5(b).  Thus,  the  discharge  characteristics  of  the 
battery  with  cross  coupling  are  evaluated  in  terms  of  the  basic 
connections. 


I 


ZL 

N2 


N3 

~r 


□_ 

N2 


N3 
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Fig.  6.  Multicell  battery  of  four  cells  in  series-parallel  connection,  (a)  With  cross 
coupling,  (b)  without  cross  coupling. 


Fig.  5.  Discharge  voltage  and  current  of  II  and  N2  in  parallel  connection  at  ambient 
temperature  of  28.5  °C.  (a)  Discharge  voltage,  (b)  discharge  current. 


4.3.2.  Without  cross  coupling 

Next,  the  experimental  results  on  the  case  of  no  cross  coupling 
are  presented  in  Table  5  and  Fig.  8.  In  Fig.  8(a),  when  II  reaches  its 
cut-off  voltage,  the  discharge  capacity  of  II  becomes  2314  mAh  and 
corresponds  to  Qs(Il).  Then,  the  discharge  capacity  of  N1  is 
restricted  by  II  and  becomes  2303  mAh.  The  discharge  capacities 
of  N2  and  N3  are  around  2700  mAh  and  smaller  than 
Qs(N2)(=Qs(N3)).  The  total  discharge  capacity  of  the  battery  is 
5009  mAh. 

By  the  same  manner  as  the  cross  coupling,  the  series-parallel 
connection  in  Fig.  6  is  aggregated  as  the  parallel  connection  of 
II  +  N1  and  N2  +  N3.  However,  any  combination  of  the  results  on 
the  basic  series  and  parallel  connections  does  not  lead  to  the 
experimentally  estimated  capacity.  Here,  the  sum  of  terminal 
voltages  of  cells  connected  in  series  is  common  for  each  string  of 
cells.  Therefore,  the  terminal  voltages  of  the  uniform  cells  N2  and 
N3  are  determined  by  the  number  Ns  of  cells  in  series  and  corre¬ 
spond  to  the  half  of  the  total  voltage  in  the  experiment.  The  ter¬ 
minal  voltage  of  N2  (or  N3)  does  not  reach  its  cut-off  voltage,  while 
the  terminal  voltage  of  II  reaches  its  cut-off  voltage.  In  this  way,  the 
discharge  capacity  of  N2  (or  N3)  does  not  reach  Qs(N2)  (or  Qs(N3)). 
Since  the  evaluation  in  the  basic  series  connection  is  based  on  the 


Table  4 

Discharge  capacities  of  four  cells  in  series— parallel  connection  with  cross 
coupling  at  ambient  temperature  of  28.6  °C. 


Discharge  capacity/mAh 

11 

2306  (23261) 

N1 

2563  (28391) 

N2 

2817  (2820f) 

N3 

2556  (2808f) 

Multicell  battery 

5121 

740 
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Fig.  7.  Discharge  voltage  and  current  of  four  cells  in  series-parallel  connection  with 
cross  coupling  at  ambient  temperature  of  28.6  °C.  (a)  Discharge  voltage,  (b)  discharge 
current. 

discharge  capacity  of  a  single  operation,  it  does  not  explain  the 
experimental  result  in  the  case  of  no  cross  coupling.  This  conclusion 
suggests  that  the  discharge  capacity  of  the  battery  depends  on  its 
configuration  and  the  terminal  voltage  of  every  cell. 

5.  Numerical  simulations 

This  section  describes  a  mathematical  model  for  the  evaluation 
of  discharge  capacity  in  the  series-parallel  connection.  Here  we 
investigate  a  dependence  of  the  battery  configurations  on  the 
discharge  capacity  of  the  general  multicell  battery  in  Fig.  9.  The 
battery  is  synthesized  with  one  ICR  cell  and  many  NCR  cells. 

5  A.  Mathematical  model 

In  Section  4.3.2,  it  is  experimentally  observed  that  the  terminal 
voltage  of  each  cell  governs  the  discharge  capacity  of  the  battery. 
Based  on  the  experimental  observation,  we  can  adopt  the  circuit- 
based  cell  model  reported  in  Refs.  [14,15]  for  analysis  of  the 


Table  5 

Discharge  capacities  of  four  cells  in  series— parallel  connection  without 
cross  coupling  at  ambient  temperature  of  29.1  °C. 


Discharge  capacity/mAh 

11 

2314  (23261) 

N1 

2303  (28391) 

N2 

2702  (28201) 

N3 

2698  (28081) 

Multicell  battery 

5009 

Fig.  8.  Discharge  voltage  and  current  of  four  cells  in  series-parallel  connection 
without  cross  coupling  at  ambient  temperature  of  29.1  °C.  (a)  Discharge  voltage,  (b) 
discharge  current. 

multicell  battery.  The  cell  model  consists  of  a  controlled  voltage 
source  and  a  constant  resistance,  denoted  as  r.  The  open-circuit 
voltage  e(t)  and  terminal  voltage  v(t)  are  represented  by  the 
following  equations: 


Ne 


Fig.  9.  Configuration  of  multicell  battery  with  nonuniform  cells.  The  battery  includes 
one  ICR  cell  and  many  NCR  cells. 
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v(t)  =  e(t)  -  ri(t), 


(2) 


1 L 


where  E0  is  the  initial  voltage,  I<  the  voltage  drop  caused  by  po¬ 
larization,  Qo  the  initial  capacity,  and  z(t)  the  discharge  current.  The 
constants  A  and  B  are  responsible  for  the  initial  voltage  drop.  The 
parameters  for  ICR  and  NCR  are  presented  in  Table  6.  They  were 
obtained  by  the  method  of  least  squares  based  on  the  experimental 
results  on  discharge  characteristics  at  constant  currents.  Here,  by 

defining  the  charge  variable  q(t)  :=  [  i(x)dx,  the  circuit-based  cell 


model  is  re-written  as 


e(<?) 

F  KQo 

0  Qo  -  <7 

d  q 

e(q)  -  v{t) 

dt 

r 

+  A  exp(-Bq), 


(3) 
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>Ns 
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The  variable  q(t)  provides  the  discharge  capacity  by  taking  t  as  a 
final  time  of  the  operation,  and  thus  the  above  circuit-based  model 
is  based  on  the  observation  in  which  the  terminal  and  open-circuit 
voltages  of  a  cell  govern  its  discharge  capacity. 

For  the  current  modeling,  we  use  the  following  two  properties 
based  on  the  experimental  results. 

PI )  The  current  of  a  battery  is  equally  distributed  to  every  par¬ 
allel  connection  with  uniform  cells. 

P2)  The  terminal  voltage  of  every  cell  is  common  in  any  series 


Fig.  10.  Classification  of  cells  in  the  multicell  battery.  All  of  the  cells  are  classified  into 
A,  B,  and  C  according  to  the  behavior  of  terminal  voltage. 

where  k  is  the  load  current.  The  matrix  formulation  of  Eqs.  (4)— (6) 
is  given  as 


connection  with  uniform  cells. 

"1  Ns  — 1 

-N5  -rA  - 

-  (Ns  —  l)rg 

Nsrc 

0 

In  the  experimental  results  on  the  series-parallel  connection  with 

0  0 

0 

1 

NP-  1 

-1 

cross  coupling,  the  current  of  the  battery  is  equally  distributed  to 

0  0 

Ns 

0 

— N5rc 

-Rl 

uniform  cells  N1  and  N3  as  shown  in  Fig.  7(b).  On  the  other  hand,  in 
the  series-parallel  connection  without  cross  coupling,  the  terminal 
voltages  of  uniform  cells  N2  and  N3  are  common  as  in  Fig.  8(a). 
These  properties  are  general  because  they  are  based  on  the  stan¬ 
dard  Kirchhoff  s  laws  for  an  electric  circuit  with  identical  cells. 

With  the  circuit-based  cell  model  and  the  above  two  properties, 
we  derive  a  mathematical  model  for  discharge  characteristics  of  the 
multicell  battery  in  Fig.  9.  Now,  all  cells  in  the  battery  are  classified 
into  three  types  of  cell  A,  B,  and  C  as  shown  in  Fig.  10.  The  classi¬ 
fication  is  based  on  the  behavior  of  terminal  voltage  of  cells. 
Because  the  current  of  every  cell  is  common  in  series  connection 
and  the  sum  of  terminal  voltages  is  common  in  the  parallel 
connection,  the  relationship  between  terminal  voltages  for  each 
string  is  described  as 


^A 

eB 

ec 

U 

h 

h 


0 

0 

0 


(7) 


eA  ~  hrA  +  (Ns  —  V(eB  —  hrs)  -  N s(ec  -  i2rc)  =  0, 


(4) 


In  a  normal  setting,  the  above  coefficient  matrix,  denoted  as 
M  =  [Mi  M2]  where  M,gR3x3,  is  column  full-rank,  and  the  sub 
matrix  M2  is  regular.  From  these,  the  following  linear  relation  of 
discharge  currents  and  open-circuit  voltages  is  obtained: 


(8) 


The  elements  Gy  of  the  above  matrix  are  presented  in  Appendix. 
Now,  based  on  Eq.  (8),  the  three  models  for  the  single  cells  A,  B,  and 
C  are  combined.  With  the  two  charge  variables  q\  and  q2  defined  as 
Jq  ij(x)dx  (j  =  1,2),  the  models  are  given  as 


~h~ 

=  -M^M, 

^A 

G11 

^12 

G13 

^A 
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eB 

= : 

G21 

G22 

G23 

eB 

k_ 

_ec_ 

G31 

G32 

G33 

_ec_ 

where  eAl  eB ,  and  ec  are  the  open-circuit  voltage  of  cell  A,  B,  and  C. 
Also,  rAl  rg,  and  rc  are  the  internal  resistance.  The  variable  i\  is  the 
current  of  cell  A  and  B,  while  i2  is  the  current  of  cell  C.  By  consid¬ 
ering  the  load  with  a  constant  resistance  RLl  the  following  equations 
are  obtained: 

z1  +  (IVp-l)i2  =  iL,  (5) 

Ns(^c  -  h rc)  -  =  0,  (6) 


Table  6 

Parameters  of  the  circuit-based  cell  model  for  two  cells  used  in  this  paper. 


eA(q 1) 


efj(Qi) 


'0,ICR 


'0,NCR 


^icrQojcr 


Qo,icr  -  Q 1 
^ncrQo,ncr 


+  AcR  exP(-#ICR<7l)> 


+  AmCR  exP(-8NCR9l) 


(9) 


(10) 


ICR 

NCR 

E0 

2.7243  V 

3.2124  V 

I< 

0.0127  V 

0.0148  V 

Qo 

2.13  Ah 

2.80  Ah 

A 

1.2006  V 

0.9475  V 

B 

0.3838  Ah”1 

0.5135  Ah1 

r 

0.1097  Q 

0.0759  Q 

Qo,ncr  -  Qi 

ec(Q2)  =  £o,NCR  —  NCR^°'N^R  +  ANcr  exP(-#NCRq2)7  (11) 

where  the  suffix  ICR  (or  NCR)  represents  the  parameters  for  the 
single  cell  ICR  (or  NCR).  Since  each  current  variable  is  represented 
as  a  superposition  of  the  three  voltage  variables  (see  Eq.  (8)),  the 
following  set  of  differential  equations  with  the  two  independent 
variables  q\  and  q2  is  obtained: 


df?l 

dt 

dq2 

dt 


=  +  GueB(q-[)  +  G-[3ec(q2), 


£21^(31)  +  G22^b(^i)  +  G23ec(q2)- 


(12) 
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These  equations  describe  the  dynamics  of  the  charges  for  the  three 
cells  A,  B,  and  C  in  a  self-consistent  manner  and  are  hence  used  as  a 
mathematical  model  for  discharge  dynamics  of  the  multicell 
battery. 

5.2.  Verification 

Here,  the  derived  model  for  the  multicell  battery  is  verified  in 
comparison  with  experiments  on  discharge  capacity  of  the  multi¬ 
cell  battery.  All  experiments  were  performed  on  the  battery  with 
one  ICR  cell  and  eleven  NCR  cells.  The  discharge  capacity  of  NCR 
was  2817  mAh  in  average,  denoted  as  Qmax,  and  was  larger  than  the 
capacity  of  ICR,  2133  mAh,  denoted  as  Qmin-  The  battery  was 
configured  as  (Ns,Np)  cells  in  series-parallel  with  the  conditions 
2  =  N$  =  6,2  =  Np  =  6,  and  N$  •  A/p  =  12. 

For  comparison,  we  calculated  the  discharge  capacity  Qof  a  cell 
in  the  series  connection  with  uniform  NCR  cells.  The  total  capacity 
of  the  multicell  battery  possibly  depends  on  the  number  of  cells, 
and  the  estimated  capacity  Q  becomes  a  quantitative  measure  of 
the  (A/s, A/p)  dependence  of  the  discharge  capacity.  The  quantity  Q  is 
simply  calculated  by  the  following  equation: 

T 

0  =  1  i2(t)dt,  (13) 

0 

where  T  is  the  final  time  of  discharge  operations.  Fig.  11  shows  the 
results  on  experiments  and  numerical  simulations  of  the  discharge 
capacity  Q  with  (A/s,A/P)  cells.  The  results  clearly  indicate  that  the 
discharge  capacity  Q  depends  on  the  number  of  cells.  In  Fig.  11, 
although  the  experimental  results  are  slightly  affected  by  the  effect  of 
cycle  life  and  a  slightly  nonuniform  characteristic  in  the  NCR  cells,  the 
numerical  and  experimental  results  are  almost  consistent.  The  derived 
model  accounts  for  the  experimental  data  on  discharge  capacity  in  the 
multicell  battery  with  nonuniform  cells.  In  this  model,  the  open- 
circuit  voltage  of  an  individual  cell  depends  on  the  amount  of 
discharge  as  well  as  discharge  capacity.  Therefore,  the  results  in  Fig.  11 
confirm  that  the  discharge  capacity  of  a  multicell  battery  is  domi¬ 
nantly  determined  by  the  open-circuit  voltages  of  individual  cells. 

In  the  experimental  results,  the  minimum  value  of  Q  is 
2416  mAh  and  the  maximum  value  is  2700  mAh.  In  the  simulations, 
the  minimum  value  of  Q  is  2409  mAh  and  the  maximum  value  is 
2701  mAh.  In  every  combination  of  Ns  and  Np,  Q  exists  in  the  range 
of  Qmin  <  Q  <  Qmax-  The  capacity  Q  depends  dominantly  on  Ns  not 
on  Np.  Also,  Q  decreases  as  Ns  increases.  The  result  on  the  discharge 
capacity  Q  <  Qmax  implies  that  the  multicell  battery  does  not 
discharge  the  full  of  its  charge.  As  the  increase  of  Ns,  this  salient 
feature  appears  clearly  and  the  multicell  battery  exhibits  a  small  Q 
as  shown  in  Fig.  11.  For  a  sufficiently  large  Ns,  from  Eq.  (4)  the  effect 
of  the  ICR  cell  A  on  the  battery  is  negligible,  and  the  battery  is 
regarded  as  the  two  NCR  cells  B  and  C  in  parallel.  Because  B  and  C 
have  the  same  characteristics,  they  will  behave  in  a  common 
manner  at  this  extreme  situation.  Here,  from  the  experimental 
result  in  Section  4.1,  the  discharge  capacity  of  the  string  containing 
A  is  determined  by  a  cell  with  the  smallest  capacity,  that  is,  A  with 
Qmin.  This  is  why  as  the  increase  of  Ns,  the  value  of  Q  decreases  and 
approaches  to  Qmin-  On  the  other  hand,  the  result  on  Q  >  Qmin 
implies  that  the  discharge  capacity  of  a  multicell  battery  is  not 
necessarily  restricted  by  a  cell  with  the  smallest  capacity.  As  dis¬ 
cussed  above,  for  the  infinite  Ns  of  cells,  the  single  ICR  cell  governs 
the  performance  of  the  battery.  On  the  other  hand,  in  a  finite  Ns,  the 
effect  of  the  ICR  cell  is  possibly  weakened  by  the  other  uniform  NCR 
cells.  Thus,  we  show  such  a  finite-number  effect  governs  the  per¬ 
formance  of  the  battery  with  nonuniform  cells.  From  the  experi¬ 
mental  result  in  Section  4.3.2,  the  finite-number  effect  is  partly 
captured  by  the  terminal  or  open-circuit  voltage  of  every  cell.  The 
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Fig.  11.  Discharge  capacity  Q.  for  the  multicell  battery  with  (NS,NP)  cells  in  the  series- 
parallel  connection,  (a)  Experimental  result,  (b)  simulation  result. 


model  derived  in  this  section  is  based  on  the  open-circuit  voltage 
and  is  therefore  successful  for  reproducing  the  finite-number  effect 
observed  in  the  experimental  results  of  Fig.  11. 

6.  Concluding  remarks 

This  paper  investigated  the  discharge  characteristics  of  a 
multicell  battery  with  nonuniform  cells.  The  analysis  reported  in 
this  paper  reveals  the  relationship  between  discharge  capacity  and 
topology  of  the  multicell  battery.  In  the  series  connection  of 
nonuniform  cells,  the  total  capacity  of  the  battery  is  governed  by  a 
cell  with  the  smallest  capacity.  In  the  parallel  connection,  the  total 
capacity  becomes  the  sum  of  individual  cells  in  parallel.  In  the  se¬ 
ries-parallel  connection  of  nonuniform  cells,  the  discharge  ca¬ 
pacity  depends  on  the  configurations  of  the  battery.  For  a  multicell 
battery  with  a  general  configuration,  the  evaluation  of  discharge 
capacity  is  enabled  by  numerical  simulations  of  the  mathematical 
model  which  is  newly  derived  in  this  paper. 

Here  we  provide  several  remarks  related  to  generalization  of  the 
experimental  results  described  above.  The  experimental  results  on 
the  series  and  parallel  connections  are  obtained  for  the  multicell 
battery  which  consists  of  less  than  two  cells.  In  the  case  of  multiple 
cells  in  the  series  connection,  the  current  of  every  cell  is  common  in 
the  same  manner  as  the  two  cells.  In  the  parallel  connection  of  mul¬ 
tiple  cells,  the  terminal  voltage  of  every  cell  is  common.  Therefore,  the 
discharge  characteristics  on  the  series  and  parallel  connections  are 
valid  in  a  multicell  battery  with  more  than  two  cells.  The  results 
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obtained  in  this  paper  lead  to  a  method  for  evaluating  discharge 
characteristics  of  multicell  batteries  with  arbitrary  configurations. 
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Appendix.  Elements  of  the  matrix  in  Eq.  (8) 
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